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Remark 3.4. The existence theoiy for diffusion approximations i ^pt entirely trivial If conduction were 
not present the question could be resolved using the well-known Kirchoff transform [S,15]. In the general 
case<I^Qre sophisticate^analysis in needed. For example the methods of [10] could be tested. 
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(^^ ' ^1 ■ ' continuous. 

A - u, 

u fc e;, 1 H\Q) n ^ ' 



Fixed point 



[ A{z{x),x)Vw-Vip = 0, V^gAVoH^)- (20) 
Jo. 

(1) z G C, C = {/-^; [Tm/n, Tmax]} closed convex. 

(3j fi v= L , L.^ is precompact. 

(^^ ' ^1 ■ ' continuous. 

A - u, 

= a{u-Ug3s). n 



□ 



Fixed point 



[ A{z{x),x)Vw-Vip = 0, V^gAVoH^)- (20) 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, £C C 

(^^ ' ^1 ■ ' continuous. 

A - u, 

= a{u-Ug3s). n 



□ 



Fixed point 



[ A{z{x),x)Vw-Vip = 0, V^gAVoH^)- (20) 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

, IS continuous. 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 

- u, 

u fc e;, 1 n c°-"(i2j. 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, 

u c e I 1 / / ^^-^^^ I 

= a{u-Ug3s). n 



Fixed point 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- (20) 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, u \s a solution. 

= a{u-Ug3s). n 



□ 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, u \s a solution. 

ue (j:nA7i(Q)nc°'"(Q). 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, u \s a solution. 

ue (j:nA7i(Q)nc°'"(Q). 

AVu -it = a{u - Ugas), a£[0,C]. 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, u \s a solution. 

ue (j:nA7i(Q)nc°'"(Q). 

AVu -it = a{u - Ug3s), ae[0,C]. 



[ A{z{x),x)VwVip = Q, V^gAVoH^)- 
Jn 

(1) z G (T, <t = {L'^;[Tmin,Tmax]} closed convex. 

(2) Maximum principles, C<L C C 

(3) (E L^, C€ is precompact. 

(4) The solution — data map, C is continuous. 
A fixed point: Cu = u, u \s a solution. 

ue (j:nA7i(Q)nc°'"(Q). 

AVu -it = a{u - Ug3s), ae[0,C]. 



Heat source 



o<f<c. vc>o, . ,r,„,„, T-,„,^ + c], 

< A(7™„, 7^3x) < A{z) < A{T^in. T^,x- C )■ 
z G e- = {Z.2 ; [7„,v,. T^,, + Q ] }, 
f 



Vw-V^= I f{z)<f. V 
Jn 



•r- 



Tmin W ^ Tmgx tin- 7fn3x)i ||^||oc)- (^4) 

Let c = c e'. 

^{Tmin- Tmax- G). < < C{Tmin- Tmax^ C)- 

C = GfC). 



Heat source 



o<f<c. iQ>o. . t;„,^ + C], 

z G e- = {Z.2 ; [7„,v,. T^,, + C ] }, 

f 



Vw-V^= I f{z)<f. V 
Jn 



•r- 



Tmin W ^ Tmgx tin- ^m3x)i ||^||oc)- (^4) 

Let c = c e'. 

^{Tmin- Tmax- G). < < C{Tmin- Tmax^ C)- 

C = GfC). 



Heat source 



0<f<C. V C > 0, , r,„,„, 7,„,^ + c ], 

z G e- = {Z.2 ; [7„,v,. 7^,, + C ] }, 



Tmin W ^ Tmgx tin- ^m3x)i ||^||oc)- (^4) 

Let C = C €. 

^{Tmin- Tmax- G). < f < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fzG[Tmin, T^ax+C], 
U ' min, ' max) ^(.^J _:: ' ^ ' m/n- ' max- )■ 

z G e- = {Z.2 ; [7„,v,. 7^,, + C ] }, 

Jn 

Tmin W ^ Tmgx tin- ^m3x)i ||^||oc)- (^4) 

Let C = C €. 

^{Tmin- Tmax- G). < < C{Tmin- Tmax^ C)- 

C = G(C). 



0<f<C. VC>0, \fz£[Tmin, Tmax+Q], 

z G e: = ; [7„,v,. 7^3x + C ] }, 

f 



Vw-V^= I f{z)<f. V 
Jn 



•r- 



Tmin W ~^max tin- ^max)i ||^||oc)- 

Let C = C €. 

^{Tmin- Tmax- G). < f < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 

z G C = {/.2 ; [7„/„, 7„,^ + C ] }, 

^m/n W ^ Tmgx tin- ^m3x)i ||^||oc)- (24) 

Let C = C €. 

^{Tmin- Tmax- G). < f < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin W ^ Tmgx tin- ^m3x)i ||^||oc)- (24) 

Let C = C €. 

^{Tmin- Tmax- G). < f < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C = ' C €. 

^{Tmin- Tmax- C). < < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C* = Ci, : c €. 

•^{Tmin- Tmax- C). < < C{Tmin- Tmax^ C)- 

C = GfC). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C = Ci, C^dt. 

•^{Tmin- Tmax- C). Q < f < C{Tmin- Tmax^ C)- 

C = G(C). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C = Ci, C^dt. 
C = GfC). 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C = Ci, C^dt. 



Heat source 



0<f<C. VC>0, \fze[Tmin, Tm^x+Q], 

< A( Tmin, 7^ax) < A{z) < A( Trr^in, T^,x, Q ), (22) 
Z G C = {/.2 ; [7™„, Trr,,^ + C ] }, 

/ >^(z)Vw-V^= / f{z)^, yip£H^{n). (23) 
Jn 

Tmin < W < Tm3x+ Ci{Q.,\[Tn,;n,Tm3x),\\f\\oo)- (24) 

Let C = Ci, C^dt. 
C* = Ci(C*). Non-existence ? 



T). L X 5, 

dtu-dh'' " ' -0. in Qt. (25) 

0, (u- 0. 
/ 'n[T^,n. 7^,,]. (27) 
A < A{z,x. t) < A. 2 G [r^/,, 7„,x]. (28) 



,7), Qr = n X S, 

dtu- f '^' ' ^0. in Qt. (25) 

0, (u- 0. 
/ 'n[T^,n. 7^,,]. (27) 
A < A{z,x. t) < A. 2 G [r^/,, 7„,x]. (28) 



,7), QT = nxS, 

^tu - div[/\(ty,x, t)Vtv] = 0, \n Qj. (25) 

0, (u- 0. 
gen''"- ''l[T^,n. 7„a,]. (27) 
A < A{z,x. t) < A. 2 G [T^i,. 7„,x]. (28) 



Parabolic model 



S = (0, 7), QT = nxS, 

^tu - div[/\(ty,x, t)Vtv] = 0, \n Qj. (25) 

(tv-^)|t=o = 0, {u-g)\3n = 0. (26) 

gen''"- 'n[T^,n. 7^,,]. (27) 

A < A{z,x. t) < A. 2 G [T^i,. 7„,x]. (28) 



Parabolic model 



S = (0, 7), QT = nxS, 

^tu - div[/\(ty,x, t)Vtv] = 0, \n Qj. (25) 

(tv-^)|t=o = 0, {u-g)\dn = 0. (26) 

g£ H\QT)n[Trr,in, T-^ax]. (27) 

A < A{z,x. t) < A. 2 G [T^i,. 7„,x]. (28) 



Parabolic model 



S = (0, 7), QT = nxS, 

^tu - div[/\(ty,x, t)Vtv] = 0, \n Qj. (25) 

(tv-^)|t=o = 0, {u-g)\dn = 0. (26) 

g£ H\QT)n[Trr,in, T-^ax]. (27) 

A < A{Z,X, t) < A, Z e [Trr^in, Tmax]. (28) 



Weak solution 



V(^G/.2(Hi), 




B. C. {u-g)e_W = {ve L^{Hl)- d 

{DtW. v3)h1(Qt) 

= [/ ^^\. 

Jn jjqt 

-.1', (y t) 



Weak solution 



{dtu, V5),2(„i)+ // A{u)Vu-Vip = 0. 



B. C. {u-g)e_W = {v e L^{Hl)] dtv G 



(29) 



{(Jtw. ^p)hhQt) 



w 



J J Qj 



','y t) 



□ 



Weak solution 



{dtu, V3)l2(h1) + / / A{u)Vu • V99 = 0. 



(29) 



B. C. {u-g)eW = {veL^{H^);dtvGL'{H-')}, (30) 



|2/u-l^ 



1) 



{(Jtw. ^p)hhQt) 



w 



J J Qj 



','y t) 



□ 



Weak solution 



V(^G/.2(Hi), 

{dtu, V3)l2(h1) + / / A{u)Vu • V99 = 0. 



B. C. {u-g)eW = {v e L^{H^y, dtv e L^{H-')} 
W-^ C([0,7-]; I.e. iu-g){x,0) = 0. 

{OfW. ip)H^QT) 

J J Qj 



(29) 

(30) 
(31) 



t) 




Weak solution 



{dtu, 99),2(„i) + jj^ A{u)Vu • V99 = 0. (29) 

B. C. {u-g)eW = {veL^iH^y,dtveL^{H-')}, (30) 

W-^ C([0,7-]; /.2(Q)), I.e. iu-g){x,0) = 0. (31) 
{dfW, '^)i.2(Hi). integration by parts? 

{OfW. ifjH^Qr) 

Jn jjqj 

-.1', (y t) 



{dtu, V3)l2(h1) + / / A{u)Vu • V99 = 0. 



Or 

B. C. {u-g)eW = {veL^{H^y,dtveL^{H-^)} 

W-^ C([0,7-]; I.e. iu-g){x,0) = 0. 

{dfW, '^)i.2(Hi). integration by parts? 



J JQj 



wdtip. 



Weak solution 



V(^G/.2(Hi), 

{dtu, V3)l2(h1) + / / A{u)Vu • V99 = 0. 



B. C. {u-g)eW = {v e L^{H^y, dtv e L^{H-')} 

W-^ C([0,7-]; I.e. iu-g){x,0) = 0. 

{dfW, '^)i.2(Hi). integration by parts? 

= [ / w^]\o - wdt^p. 

Jo. JjQ-r 
1 /■*+'' 

Wh{x,t) = w{x,s)ds. 



(29) 

(30) 
(31) 

(32) 
(33) 



Linearized map 



e = {0 ); Tmin < (p{x, t) < Tmax^ a. e. in Qt}- 



{w-g)(^W. (m/-^,, (36) 
Cz = w. C€ C e*? 




Linearized map 



(t = {(^e L^{Qt)\ Trr^in < H^, t) < Trr,,^, a. e. in Qj] . (34) 

Vy - J, 



(flfw, v5)^2(S;Hoi(n)) + jj^ A{z)Vw 



{w-g)(^W. (m/-^,, (36) 
Cz = w. C€ c e*? 



Linearized map 



€ = {(b£ L^{Qt): T^in < Hx, t) < T^ax, a. e. in Qj}. (34) 
is a closed convex set in L^{Qt)- 



iw-g)eW. {w- - - (36) 
£z = w. C€ C e*? 



Linearized map 



= L^{Qt); T^in < H^, t) < T^sx, a. e. in Qj}. (34) 
is a closed convex set in L^{Qt)- 

iw-g)eW. (m/-^,, (36) 
£z = w. C€ C e*? 




Linearized map 



= L^{Qt); T^in < H^, t) < T^sx, a. e. in Qj}. (34) 
is a closed convex set in L^{Qt)- 



(5tw, ^)^2(S;Hi(n)) + / / A{z)VwV^ = 0. (35) 







T 



iw-g)eW. {w- ... (36) 
£z = w. C€ C e*? 




Linearized map 



= L^{Qt); T^in < H^, t) < T^sx, a. e. in Qj}. (34) 

is a closed convex set in L^{Qt)- 
V(pG/.2(S; Ayi(Q)), 



(5tw, ^)^2(S;Hi(n)) + / / A{z)VwV^ = 0. (35) 







T 



{w-g)GW, {w-g){x,0) = 0. 

Cz = w. 



(36) 




Linearized map 



= L^{Qt); T^in < H^, t) < T^sx, a. e. in Qj}. (34) 

is a closed convex set in L^{Qt)- 
V(pG/.2(S; Ayi(Q)), 



(5tw, ^)^2(S;Hi(n)) + / / A{z)VwV^ = 0. (35) 







T 



{w-g)GW, {w-g){x,0) = 0. 

Cz = w. 



(36) 




Linearized map 



= L^{Qt); T^in < H^, t) < T^sx, a. e. in Qj}. (34) 
is a closed convex set in L^{Qt)- 



(5tw, ^)^2(S;Hi(n)) + / / A{z)VwV^ = 0. (35) 







T 



{w-g)GW, {w-g){x,0) = 0. (36) 
Cz =w. etc a 



Maximum principles (1) 



y{w-g)eW. 3{v/} c C~(Hi), --u,||vv- 

{dtivi + g). [Vi+g I /„^x) + >/.2(Hi) 







\2 



[ + W+(x. 7-) _ ^ 
Jn ^ 



Maximum principles (1) 



y{w-g)eVV. 3{v/} c C-(Hi), --u,||vv- 

{dtivi + g). [Vi+g I /„^x) + >/.2(Hi) 







\2 



[ + W+(x. 7-) _ ^ 
Jn ^ 



Maximum principles (1) 



W = {v e L^H^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 
y{w-g)eVV, 3{vi} c C°^{H^), --ViWw^O, 

{dtivi + g). {Vi+g I /„^x) + >/.2(Hgi) 



{0. 





rr 



"""" '0+ 



0. 



(38) 



Maximum principles (1) 



W = {v e L^iH^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 



y{w-g)eW, J{u/}cC-(Hi) 



- Vi \W 



{Otivi + g). [Vi +g I /„^x) + >/.2(Hgi) 



{0. 





rr 



"""" '0+ 



0. 



(38) 



Maximum principles (1) 



W = {v e L^iH^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 



y{w-g)eW, 3{vi}cC^{H^) 



{0. 





rr 



- vi\\w 0, 

I /iwx) + >/.2(Hgi) 

"""" '0+ 



0. 



(38) 



Maximum principles (1) 



W = {v e L^iH^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 
y{w-g)eW, 3{u/} C C-(AVi), \\w-g- viWw^O, 



{0. 





rr 



"""" '0+ 



0. 



(38) 



Maximum principles (1) 



W = {v e L^iH^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 
y{w-g)eW, 3{u/} C C-(AVi), \\w-g- viWw^O, 

{dt{vi + g), {Vi + g- "rmax)+) z.2(„l) 

= ^ {dt{vi + g), {vi + g- 

dt{vi + g) - {vi + g - Tm3x)+ (37) 



■J -J 



0. 



(38) 



Maximum principles (1) 



W = {v e L^iH^y, dtv e L^{H-^)} = C-([0, 7]; AVi(n)). 
y{w-g)eW, 3{u/} C C-(AVi), \\w-g- viWw^O, 

{dt{vi + g), {Vi + g- "rmax)+) z.2(„l) 

= ^ {dt{vi + g), {vi + g- 

dt{vi + g) - {vi + g - Tm3x)+ (37) 



dt 



{vi + g- Tmax)\{x, i) 



{v;+g - 7max) + (x, T) 



0. 



(38) 



Maximum principles (2) 



{dtW, {W - r^ax) + )L2(Hl) = ^ > 0- (39) 



J Jqj 

:nax)- 



> 



J J Qt 

r r 



> . .. , , I {w- T„ax)i. (40) 

Js Jn 



Maximum principles (2) 



{dtW, {W - r^ax) + )L2(Hl) = ^ > 0- (39) 



A{z)VwV{w- 7^,^)+ 

A{z)V{w - 7max)+ • V(m/ - Trr^ax)- 



T 

> C(Q)A [ [{w- Tn.,x)\. (40) 
JS Jo. 



Maximum principles (3) 



let ip = {w - Tmax) ■ , ,). 

= (t'tlV. (w - T^ax) + )/.2(s.H„i(n)) 

> + ' - ' ' - T^,,)l > 0. (41) 

IKIIw = IKIIl2(Ho1) + I' ' ' " ,4-1) ^ (42) 

W is precompact. 



Maximum principles (3) 



Let ip = {w- Tmax)+ G /-^(^o), 

> + ' ' - T^,,)l > 0. (41) 

IKIIw = IKIIl2(Ho1) + I' ' ' " ,4-1) ^ (42) 

W is precompact. 



Maximum principles (3) 



Let ip = {w- Tmax)+ G L'^iH^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ A{z)Vw-V{w- T^,x)+ 
J Jqj 

> + C(Q)A I I [w- T^,x)l > 0. 
Js Jn 



W < Trr 



Kllvv 



is precompact. 



\w 



(41) 



(42) 



Maximum principles (3) 



Let ip = {w- Tmax)+ G L'^iH^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ A{z)Vw-V{w- T^,x)+ 
J Jqj 

> + C(Q)A I I [w- T^,x)l > 0. 
Js Jn 



w < T„ 



llWllW 



is precompact. 



2 _ ||..,||2 



\W\ 



(41) 



(42) 



Maximum principles (3) 



Let ip = {w- Tmax)+ G L'^iH^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ A{z)Vw-V{w- T^,x)+ 
J Jqj 

> + C(Q)A I I [w- T^,x)l > 0. 
Js Jn 



(41) 



W < Tr, 



+ 



(42) 



is precompact. 



Maximum principles (3) 



Let ^ = [w- T^,^)+ G L^{H^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ [[ A{z)VwV{w- Trr^^xY 



w < T„ 



> + C(Q)A [w- T^,,)\ > 0. 
Js Jn 

|2 , HQ ...ii2 



Kllvv 



1/1/ 



is precompact. 



(41) 



(42) 



Maximum principles (3) 



Let ^ = [w- T^,^)+ G L^{H^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ [[ A{z)VwV{w- Trr^^xY 



w < T„ 



> + C(Q)A {w- 7„,,)| > 0. 
Js Jn 

|2 , HQ ...ii2 



Kllvv 



1/1/ 



is precompact. 



(41) 



(42) 



Maximum principles (3) 



Let ip = {w- Tmax)+ G L'^iH^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ A{z)Vw-V{w- T^,x)+ 
J Jqj 

> + C(Q)A I I [w- T^,x)l > 0. 
Js Jn 



|2 , HQ ...ii2 



W (E /-^(Qr). -Ce: is precompact. 



(41) 



(42) 




Maximum principles (3) 



Let ip = {w- Tmax)+ G L'^iH^). 

= {dtW, {w - Trr^3x)+)L2{S■H^{n)) 

+ A{z)Vw-V{w- T^,x)+ 
J Jqj 

> + C(Q)A I I [w- T^,x)l > 0. (41) 
Js Jn 

IKIIvV = IKIIl2(h1) + ll^t'^llL2(H-i) ^ ("^2) 

W (E L^{Qt)^ >Ce: is precompact. 
A fixed point \^ C is continuous. 



C is continuous. 



,4(z)Vm/ ■ 



0. 



0. (44) 



0, 



(St wo- V^)z.2(Hi) 



f^t M// ^ ' 

Vw; ^ Vm/q, 



Vip = 0. 



'0 • V 



(45) 



W/ — VI/ 1 2 



L is continuous. 



V(^G z/, zee:, 



0, 



(St wo- V^)z.2(Hi) 



f f A{z)Vw ■ 
dt Wj dt i 



0. 



0. (44) 



'0 • Vv? = 0. 



(45) 



W/ — VI/ 1 2 



L is continuous. 



V(^G z/, zee:, 

{dtWi,ip)L2^H^^) + JJ^ A{zi)VwrVip = 0, (43) 



WAizi 



{dtw, ^)l2{h^) +11 A{z)Vw ■ Vip 
■ 0, 



(44) 



dt Wj dt i 
Vwi Vm/q, 



/.2. 



Vip = 0. 



'0 • V 



(45) 



W/ — VI/ 1 2 



L is continuous. 



V(^G z/, zee:, 

{dtWi,ip)L2^H^^) + JJ^ A{zi)VwrVip = 0, (43) 



{wi} is bounded in W. 

\\A{zi) - , • 0, Vwi Vm/o, V 



(44) 



/.2. 



7 



Vip = 0. 



'0 • V 



(45) 



W/ — VI/ 1 2 



L is continuous. 



V(^G z/, zee:, 

{dtWi,ip)L2^H^^) + JJ^ A{zi)VwrVip = 0, (43) 



((?tM/, (/')z.2(Hi) + J J A{z)\/w-\/ip = 0. 
{m//} is bounded in W. f?ti/i// ^ dtWQ, weakly in 
\\A{zi) - A{z)\\2 0, Vwi Vwq, weakly in L^. 



(44) 



7 



Vip = 0. 



■ ■ V 



(45) 



W/ — VI/ 1 2 



L is continuous. 



V(^G z/, zee:, 

{dtWi,ip)L2^H^^) + JJ^ A{zi)VwrVip = 0, (43) 



((?tM/, (/')z.2(Hi) + J J A{z)\/w-\/ip = 0. 
{m//} is bounded in W. f?ti/i// ^ dtWQ, weakly in 
\\A{zi) - A{z)\\2 0, Vwi Vwq, weakly in L^. 



(44) 



(9tM/o, '^)z.2(Hi) + //^ /\(z)Vwo • = 0. (45) 



L is continuous. 



V(^G z/, zee:, 

{dtWi,ip)L2^H^^) + JJ^ A{zi)VwrVip = 0, (43) 



((?tM/, (/')z.2(Hi) + J J A{z)\/w-\/ip = 0. 
{m//} is bounded in W. f?ti/i// ^ dtWQ, weakly in 
\\A{zi) - A{z)\\2 0, Vwi Vwq, weakly in L^. 



(44) 



(9tM/o, '^)z.2(Hi) + //^ /\(z)Vwo • = 0. (45) 



Wq = W. W (£ L^{Qt), \\wi — w\\2 — )• 0. 



Integral 



dtw 



1 n 



. -ho 



z(x,s)ds)VH/] = 0. (46) 



(^tw, v?>i.2(/-, z(s)di 

" J JQt "0 Jt-ho 

Sir ' - 



Integral 



If t7|[_ftg 0] is known, 



1 rf 

dtw - z{x.s)ds)Vw] = 0. (46) 

. -ho 



(^tw, v?>i.2(/-, z(s)df, 

J JQt "0 Jt-ho 



Sir 



Integral 



If t7|[_ftg 0] is known, 

1 r 



dtw-d\v[A{^ [ z(x,s)ds)VH = 0. (46) 

"0 Jt-ho 



J JQt "0 Jt-ho 

Sir ' - 



Integral 



If t7|[_ftg 0] is known, 

1 f 

dtw-d\v[A{— z(x,s)ds)VH = 0. (46) 

"0 Jt-ho 

{dtw,ip),2^Hh+ ff I z{s)ds)Vw-V^ = 0. (47) 

J jQr "0 Jt-ho 

Sir ' 



If t7|[_ftg 0] is known, 

1 f 

dtw-d\v[A{— z(x,s)ds)VH = 0. (46) 

"0 Jt-ho 

{dtw,ip)i_2(HU+ ff [ z{s)ds)Vw-V^ = 0. (47) 

°' J jQr "0 Jt-ho 

Similarly, 



Elastic wave 




l_oc 7 



Elastic wave 



Umax is independent of A, A. 



Elastic wave 




Umax is independent of A, A. 

+ //^/S'(.x,0gg^0. (48) 



Elastic wave 




Umax is independent of A, A. 

+ //^/S'(.x,0gg^0. (48) 

L°° ? 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 



Su, , . ^ , (z,x. t) e /- , z, e.g. f e /."^. 

(2) n = l, C(Ayi)- 

(3) /\(x.t) e [A, A], zee-,. (Or); [-Q. Q]}. 
£0" ^ o\ 

j precompact. 



Wave 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Su, , , , (z,x. t) el'- y. e.g. f e /."^. 

(2) n = l, C(Ayi)-- 

(3) /\(x.t) e [A, A], zee-,. (Or); [-Q. Q]}. 

£0" 

j precompact. 



Wave 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(■ , 

(3) A{x.t) G [A, A], zG e - (Or); [-Q. G]}. 
£0" ^ o\ 

j precompact. 



Wave 



{w'\^)i_2^Hl) + jj^ A{x,t)VwVip = jj^ f{z,x,t)ip. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 
(2) n = 1, C(Ayi) ^ C%Qr). 

. . . . L . / j ; [-Ci, Ci]}. 

£0" ('■. 

j precompact. 

is ^npti'- 




Wave ... 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

'^llc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(2) n = 1, C(Ayi) ^ CO(Qr). 

(3) /\(x,t)G[A,A], zG£ = {/.2(Qr); [-Q, Q]}. 

j precompact. 



is ^npti'- 




Wave ... 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(2) n = 1, C(Ayi) ^ CO(Qr). 

(3) /\(x,t)G[A,A], zG£ = {/.2(Qr); [-Q, Q]}. 

H^(Or} j precompact. 



is ^npti'- 




Wave ... 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(2) n = 1, C(Ayi) ^ CO(Qr). 

(3) /\(x,t)G[A,A], zG£ = {/.2(Qr); [-Q, Q]}. 

H^{Qt) (e L^{Qt)- - precompact. 



is ^npti'- 




Wave ... 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(2) n = 1, C(Ayi) ^ CO(Qr). 

(3) /\(x,t)G[A,A], zG£ = {/.2(Qr); [-Q, Q]}. 

H'^{Qt) (e /-^(Qt). -Ce: is precompact. 



is '-nnti'- 




Wave ... 



(^", V.),2(Hi) + jj^ A{x,t)VwV^ = jj^ f{z,x,t)^. (49) 

lkllc(Hi) + lk'llz--(i.2) + lk"llz.2(H-i) < Ci(||/'||2 + ||(?!>o||h1 + II</'i|Iz_2)- 

(50) 

Suppose (1) f{z,x, t) G for any z, e.g. f G /-°^. 

(2) n = 1, C(Ayi) ^ CO(Qr). 

(3) /\(x,t)G[A,A], zG£ = {/.2(Qr); [-Q, Q]}. 

Hi((?r) (e L^iQr)- C€ is precompact. 
A fixed point if C is continuous. 



Part 2: Math algorithms 



Qt 



J J Qi 



(ptu. 



a:(u;+i - Ugss)^ = II f{u;)ip. 
(0,t)xr JjQt 



A{ui. , /,). 

A{ui+,)^A{ui) + A'{ui){ui^i~Ui), ^' 

0: 



Part 2: Math algorithms 



Qt = nx (0, f). 



+ J 1^ A{u,)Vui+i-V^ 



(o,t)xr J./Qt 



A{ui. , /,). 

/\(u/+i)«/\(£y,) + /\'(u/)(u/+i-u,-), ^' 

0: 



Part 2: Math algorithms 



Qt = Q. X (0, t). If t is small enough, 

+ "(tv;+i - Ug^s)^ = / / f{ui)ip. 

././(o,t)xr J./Qt 

A{ui+,) « /\(£y/) + /\'(u/)(u/+i - ui), ^ - " - H. 

(I 



Part 2: Math algorithms 



Qt = Q. X (0, t). If t is small enough, 



// A{ui)Vui+i-V^ 

J jQt 



+ a(tv/+i - Ugas)ip = / / f{ui)ip. 

JJ{o,t)xr JJQt 

A{ui+i) K A{ui). 

A{ui+,) « A{ui) + /\'(u/)(u/+i - ui), 

(I 



(51) 



Part 2: Math algorithms 



Qt = Q. X (0, t). If t is small enough, 



// A{ui)Vui+i-V^ 

J jQt 



+ a(tv/+i - Ugas)ip = / / f{ui)ip. 

JJ{o,t)xr JJQt 

A{ui+i) « A{ui). 

A{ui+,) ^ A{ui) + /\'(u/)(u/+i - ui), 

(I 



(51) 




Part 2: Math algorithms 



Qt = Q. X (0, t). If t is small enough, 

+ a(tv/+i - Ugas)ip = / / f{ui)ip. (51) 

JJ{o,t)xr JJQt 

A{ui+i) ^ A{ui). 

/4(u/+i) ?a A{ui) + /4'(t7;)(tv/4_i — u/), Newton method. 

0: 



Part 2: Math algorithms 



Qt = Q. X (0, t). If t is small enough, 



// A{ui)Vui+i-V^ 

J jQt 



+ a(tv/+i - Ugas)ip = / / f{ui)ip. 

JJ{o,t)xr JJQt 

A{ui+i) ^ A{ui). 

/4(u/+i) ?a A{ui) + /4'(t7;)(tv/4_i — Uj), Newton method. 
_ Ro^hg method. 



(51) 



m 



Part 3: SOTS — homogenization 



u,{x) = uo{x) + eui{x. -) + s^U2{x. -) + ... 



,. x) + sAi{uo, ui. X, -) + s^A2 + 



a°j{uo.x) 



J 



■Vif = 0. 

;>o.^,y)]dy. 



(53) 
(54) 
(55) 



m 



Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 

u,(x)="^ ^ '-sui{x.-) + e^U2{x.^) + ... (53) 

,. x) + sAi{uo, ui. X, ^) + s^A2 + ... (54) 

• Vv? = 0. (55) 

f 

3°{uo.x) = ^ \uo.x,y)]dy. 



m 



Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 

tve(x) = tvo(x) + eui{x, ^) + e^U2{x, ^) + ... (53) 

,. x) + eAi{uo, ui. X. ^) + s^A2 + ... (54) 

• Vv? = 0. (55) 

f 

3°{uo.x) = ^ \uo.x,y)]dy. 



m 



Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 

Ue{x) = tvo(x) + eui{x, ^) + e^U2{x, ^) + ... (53) 

A{ue) = A{uo,x) + £Ai{uo, tvi,x, ^) + e^A2 + ... (54) 

• Vv? = 0. (55) 

f 

3°{uo.x) = ^ \uo.x,y)]dy. 



Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 

Ue{x) = tvo(x) + eui{x, ^) + e^U2{x, ^) + ... (53) 

A{ue) = A{uQ,x) + £Ai{uo, tvi,x, ^) + e^A2 + ... (54) 

/ Ao{uo,x]Vuo-Vip = 0. (55) 
Jn 

3°{uo.x) = ^ \uo.x,y)]dy. 




Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 
Ue{x) = tvo(x) + eui{x, ^) + e^U2{x, ^) + ... (53) 
A{ue) = A{uQ,x) + £Ai{uo, tvi,x, ^) + e^A2 + ... (54) 

/ Ao{uo,x]Vuo-Vip = 0. (55) 
Jn 

f d 
a?(tvo,x)= / [aij{uo) + aii{uo)^Nj{uo,x,y)]dy. (56) 

Jy (jyi 



Part 3: SOTS — homogenization 



-d\y[A{ue,x,-)Vue] = 0. (52) 

Ue{x) = tvo(x) + eui{x, ^) + e^U2{x, ^) + ... (53) 

A{ue) = A{uQ,x) + £Ai{uo, tvi,x, ^) + e^A2 + ... (54) 

/ Ao{uo,x]Vuo-Vip = 0. (55) 
Jn 

f d 

a-!(tvo,x)= / [aij{uo) + aii{uo)^Nj{uo,x,y)]dy. (56) 

Jy (jyi 

Aiuo)V{Nj + yj) • = 0, V(/7 E W^.^iY). (57) 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

1 la, A < Ao{z) < Ao. 

If /4(z) is cu, .1 ^.ous, '^o(^) is ^ iivj . .J iL i.jous. 

3uo. s.t. 7uQ-V^ = 0. (58) 

■hi 

II - "o|lc''{n) ~^ 0' 

2 

Fr - ^ A. - " ^))/.,(x)/..y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

1 la, A < Ao{z) < Ao. 

If /4(z) is cu, .1 ^.ous, '^o(^) is ^ iivj . .J iL i.jous. 

3uo. s.t. 7uQ-V^ = 0. (58) 

■hi 

II - "o|lc''{n) ~^ 0' 

2 

Fr - ^ A. - " ^))/.,(x)/..y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Aq{z) < Aq. 

If /4(z) is cu, .i.i^.ous, '^o(^) is ^ iivj . .J iL i.jous. 

3uo. s.t. 7uQ-V^ = 0. (58) 

■hi 

II - "o|lc''{n) ~^ 0' 

2 

Fr - ^ A. - " ^))/.,(x)/..y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao[z) < Aq. 

If /4(z) is cu, .i.i^.ous, '^o(^) is ^11- . .J iL i.jous. 

3uo. s.t. 7uQ-V^ = 0. (58) 

■hi 

II - "o|lc''{n) ~^ 0' 

2 

Fr - ^ A. - " ^))/.,(x)/..y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao{z) < Aq. 

If /4(z) is continuous, , ,u^..y . ..lo., . u u i.jous. 



3uo. s.t. 7uQ-V^ = 0. (58) 



\ue - "ol!c''{n) 



0, 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao{z) < Aq. 

If A{z) is continuous, '^o(^) is also continuous. 

3uo. s. t. 7uQ ■ Vi/J = 0. 

■hi 

II - "o|lc''{n) ~^ 0' 

2 

Fi " . ' ^))/.,(x)/..y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao{z) < Aq. 

If A{z) is continuous, '^o(^) is also continuous. 



3uo, s.t. 



/ Ao{uQ,x)Vuo-Vip = 0. 
Jn 



Aiu ^ ' (>^)/V-y)- 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao{z) < Aq. 

If A{z) is continuous, '^o(^) is also continuous. 

3uQ, s.t. / /4o(ivo,^)Vuo • V(/? = 0. 
Jn 

II - "ollco(n) ~^ 0' 

/4(tV£)VtV£ ^ Ao{uo)Vuo, weakly in L^. 

Fi - ' Ai^ - ■ W)),v.(x),v.y). 



Homogenization (1) 



lfze[7™v,, 7„3x], X<A{z)<A. 

The div-curl Lemma, A < Ao{z) < Aq. 

If A{z) is continuous, '^o(^) is also continuous. 

3uo, s.t. / Ao{uQ,x)Vuo-Vip = 0. (58) 
Jn 

II - "ollco(n) ~^ 0' 

/4(ue)VtVe ^ Ao{uo)Vuo, weakly in L^. 

Fusco, Moscariello, /4,v = /A((m/(x)),v, (x),v, y). 



Homogenization (2) 



Ve > 0, 3u-r, s.t. 



■ V93 = 0. II Uj.^ — t^ollc" ^ 0- 



/ A{uoA)Vuo,s-Vip = 0. [ Ao{u, , 

\\u0.e - W0II/.2 -> 0, A{uq. -)Vuo.t- - 

||u, - uollc" ^ 0. f ) - ^("0, f )lk- ^ 0, 

II — Uo.ellfyl — > 0. II Uj. — Uo.t-||i^2 — > 0. Uq = Wq. 



Homogenization (2) 



Ve > 0, 3 Ue, s.t. 

I A{u,,-)Vus-V(p = 0. -uolico ^0. (59) 

Jn e 

/ A{uoA)Vuo,s-Vip = 0. I Ao{u, , 

\\u0.e - W0II/.2 -> 0, A{uq. -)Vuo.t- - 

- uollc" ^ 0. f ) - ^("0, f )lk- ^ 0, 

II — Uo.ellfyl — > 0. II Uj. — Uo.t-||i^2 — > 0. Uq = Wq. 



Homogenization (2) 



Ve > 0, 3 Ue, s.t. 

I A{u,,-)Vus-V(p = 0. -uolico ^0. (59) 

Jn e 

/ /4(uo,-)Vuo.£ • V(^ = 0, / Ao{uo)Vwo-Vip = 0. (60) 
Jn ^ Jn 

\\uo.e - M/0II/.2 -> 0, /^(tVo. -)Vuo.t- - 

- uollc" ^ 0. f ) - ^("0, f )lk- ^ 0, 

II — Uo.ellfyl — > 0. II Uj. — Uo.t-||i^2 — > 0. Uq = Wq. 



Homogenization (2) 



Ve > 0, 3 Ue, s.t. 

I A{u,,-)Vus-V(p = 0. -uolico ^0. (59) 

Jn e 

/ /4(uo,-)Vuo.£ • V(^ = 0, / Ao{uo)Vwo-Vip = 0. (60) 
Jn ^ Jn 

\\uo,£ - wo\\l2 ^ 0, A{uo,-)\/uo^£ ^ Ao{uo)\/wo. (61) 

- tvollc" ^ 0, \\A{u,. ^j . .v"o, f )||l- 0, 

II — Uo.ellfyl — > 0. II Uj. — Uo.t-||i^2 — > 0. Uq = Wq. 



Homogenization (2) 
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